Centrifugal terms in the WKB approximation and semiclassical quantization of 

hydrogen 
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A systematic semiclassical expansion of the hydrogen problem about the classical Kepler problem 
is shown to yield remarkably accurate results. Ad hoc changes of the centrifugal term, such as the 
standard Langer modification where the factor l{l + 1) is replaced by {I + 1/2)^, are avoided. The 
semiclassical energy levels are shown to be exact to first order in h with all higher order contributions 
vanishing. The wave functions and dipole matrix elements are also discussed. 
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While the solution of the hydrogen problem was one of 
the early successes of quantum mechanics, it failed to be 
a showpiece for the WKB approximation which proved 
to do rather poorly. Usually, this is attributed to the 
singularity of the Coulomb potential at r = where r is 
the distance between proton and electron. Clearly, near 
the origin the WKB expansion cannot be justified even 
in the semiclassical limit. Langer |l] has shown that the 
correct behavior for r — > can be enforced to the WKB 
wave function if the centrifugal term 



Vc{r) 



Tt?1{1 



1) 



2mr^ 



in the radial Schrodinger equation is replaced by 



(1) 



(2) 



Quite remarkably, with the Langer modification (LM) 
of the interaction potential, the WKB approximation 
gives exact energy eigenvalues for the hydrogen problem 
already to lowest order. As a consequence, the LM is 
now seen as a standard ingredient of WKB theory for 
the hydrogen problem and related systems with radial 
symmetry, such as the radial harmonic oscillator or the 
Morse potential in three dimensions. For recent applica- 
tions and extensions we refer to the work by Yi et al. ||] 
and by Morehead |^ . 

In the last years some attempts have been made to 
avoid the Langer modification. For the exactly solvable 
hydrogen problem semiclassical theories based on nonlin- 
ear transformations jj] or supersymmetry Q are power- 
ful alternatives to conventional WKB methods. However, 
these approaches lead to exact results only for the strict 
1 /r potential and do not constitute a general replacement 
for the standard semiclassical expansion. Within the con- 
ventional approach, Friedrich and Trost Q| have avoided 
the LM, introducing instead an additional phase of the 
WKB wave function which is then optimized. For the 
repulsive potential their method give results that 
are superior to those derived from conventional WKB 



with LM. However, also their approach maintains that 
for Coulomb type problems the textbook WKB expan- 
sion needs to be modified. In this Letter we challenge 
this common believe. 

We start from the obvious observation that in the clas- 
sical limit the hydrogen problem should reduce to the 
Kepler problem. The form of the classical orbits depends 
on the energy E and the angular momentum L. Hence, 
the leading order WKB radial wave function should also 
be calculated for given E and L = til. This implies that 
within the WKB expansion the centrifugal potential term 
(|l|) should be decomposed as 



Vc{r) = 



L 



2mr- 
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where the first term is the classical centrifugal term while 
the second term is a quantum correction. Since the WKB 
expansion proceeds in powers of h, this latter term has 
to be treated as a perturbation and expanded accord- 
ingly. Remarkably, the consequences of such a strictly 
systematic expansion in powers of h seem not to have 
been investigated previously. 

We demonstrate that a systematic expansion (SE) 
about the Kepler problem yields WKB wave functions 
that are as accurate as for other potential problems de- 
spite the singularity at r = 0. Notably, the semiclassi- 
cal energy eigenvalues for the hydrogen problem become 
exact to first order in h with all higher order correc- 
tions vanishing, while for the problem with LM the exact 
semiclassical eigenvalues obtained in lowest order become 
worse when higher order corrections are evaluated |^]. 

We start from the radial Schrodinger equation for the 
hydrogen atom 



2m dr"^ 



e 
r 



Vc{r) U-W = E^{r) 



(4) 



with Vc{r) given by Eq. (||). Using the conventional 
WKB ansatz for the wave function 



1 



^'(r) — exp 
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k=0 



and expanding in powers of h, we obtain for the quanti- 
ties 



yk{r,E,L)^dSkir,E,L)/dr 
the recursive set of equations 
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and where p{r,E,L) — yo{r,E,L) is the classical mo- 
mentum. Further, the prime denotes differentiation 
with respect to r. These equations yield two functions 
y^^\r,E^L) depending on the choice of the sign of the 
momentum p{r,E,L), and the wave function is a linear 
combination of the form 




^{r,E,L) = c^'^^exp 

CT = ± 

where 

oo 

y(r,i?,L) = ^(-z;i)'=2;fe(r,£;,L) 



dry'^''\r,E,L)\ (12) 
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The momentum p{r, E, L) has a branch cut which is cho- 
sen conveniently between the classical turning points 



?'i,2 = a (1 T e) : 



(14) 



where a is the big axis and e the eccentricity of the ellipse 
in the Kepler problem. Dunham ||] has shown that by 
choosing the initial point of integration ro on the left side 
of the two classical turning points and a contour avoid- 
ing the turning points as indicated in Fig. la, the wave 
fmiction becomes 



* (r, E, L) 
with 



elsewhere 
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^'^^\r,E,L)=exp(-J dr y^^\r, E, L) ] . (16) 
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FIG. 1. a) The complex r plane with the classical turning 
points ri/2- Connecting points of the classical allowed and 
forbidden regions one has to avoid the turning points by in- 
tegrating along the circles, b) Deformation of the integration 
contour in the complex plane. 

Since we search for a unique solution, we have to re- 
quire that the wave function is independent of whether 
one integrates above or below the branch cut. This leads 
to the condition 



— ^ dry (r, e, L) = 27ri 



[Ur + 1) ; 



(17) 



where rir is a positive integer and the integration contour 
encircles the branch cut. Using this equation one gets a 
quantization of the energy which is related to the Bohr- 
Sommerfeld rule. To evaluate the contour integrals, we 
use a technique due to Sommerfeld which exploits the 
fact that the yk {r, E, L) have only poles on the positive 
real axis. By this assumption we find As indicated in Fig. 
lb one has to calculate integrals along the contours C2 
and C3 instead of encircling the branch cut. To order h 
the integrals are readily evaluated yielding 



2Eh^ 



1, 



which gives the exact energy eigenvalues for the bound 
states of the hydrogen atom 



En 
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-me /2h n 
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with the principal quantum number n = Ur + 1 + I. Cor- 
rections of higher order in h coming from the contour 
integrals over the functions ykj k > 2 vanish exactly. To 
show this we first investigate the analytical structure of 
yo and yi at the origin. We find 



— > 

Re[r] 
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,E,L) = iLr-^ + 0{r°) 
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while the power series expansion of yi begins with a hn- 
ear term. Consequently the expansion of y[ starts with 
a constant term. Now, using 



2/2 



yl + y'l 
22/0 



(20) 



one immediately sees that the expansion of 2/2 begins with 
a linear term and therefore the residue of j/2 at the ori- 
gin is zero. Since the recurrence relations ^) and ( |l0| ) 
contain j/o only in the denominator, it is easy to show 
by induction that the Taylor series of all yk with k > 2 
start with linear or higher order terms. This implies that 
the integrals along the contour C'2 vanish for all yk with 
k > 2. In an analogous way one can treat the integrals 
along the contour C3 by replacing r by 1/u and remem- 
bering the additional factor — 1/w^ originating from the 
transformation of the integration measure. One finds 
that the integrals along the contour C3 also vanish for 
all yk with k > 2. Therefore the semiclassical energy 
quantization (^) is exact to all orders in h, while in the 
WKB approximation with Langer modification higher or- 
der terms destroy the exactness of the energy eigenvalues. 

Next we consider the wave functions. Disregarding 
quadratic and higher powers in h in (|l^) and (p^, we 
arrive at an expression for the lowest order WKB wave 
functions for r on the positive real axis of the form 



¥-\r,E,L) + ^^+\r,E,L) (21) 



with 



(22) 



where the additional phase ip{r, E, L) arises from the part 
of the centrifugal term in (0) that is linear in Ti. In fact. 



^{r,E,L) = -dS„{r,E,L)/dL 



(23) 



is just the phase of the classical trajectory in the plane of 
motion of the Kepler problem in terms of which eq. (|^) 
can be written as 
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A representation of the WKB wave function as the real 
part of the superposition of incoming and outgoing waves 
as in Eq. (|2l| ) was introduced previously by More and 
Warren for the standard approach with LM. Since 
the undesirable growing part of the wave function has 
a purely imaginary coefficient, it is removed when the 
real part is taken. The normalization c{E, L) of the wave 
function is obtained from 



/ri 

which gives 



r2 1 

dr ^ -(j) dr = 1. (25) 



2m dE 
nh^ dn ' 



(26) 



More and Warren refer to the omission of the terms 
+ v]/(-)vl/(-) in the normalization integral as 
" restricted interference approximation" . Finally, we get 
for the WKB wave function of the hydrogen atom in the 
classical accessible region between the two turning points 



■^(r,E,L) = ^^cos 

Vp 
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We now compare the WKB wave functions with the ex- 
act ones. A typical feature of WKB wave functions is 
the divergence at the classical turning points. As can be 
seen from Fig. 2, this behavior is qualitatively the same 
for the Langer modified expansion (LM), our systematic 
?i-expansion WKB(SE), and poor man's WKB(PM) ob- 
tained when the full centrifugal term (|^) is retained in 
the lowest order equation. While the WKB(PM) wave 
function for the ground state does indeed poorly, the 
main difference between the WKB(LM) and WKB(SE) 
wave functions comes from the fact that the distance be- 
tween the turning points of the Langer modified wave 
functions is smaller. This is just a consequence of the 
shift of the turning points due to the LM. Therefore, be- 
tween the turning points, our wave functions give a bet- 
ter approximation to the exact ones. For the s-states, the 
Langer modified wave functions are constructed to vanish 
at r = and they have a divergence near the origin since 
the left turning point is moved away from r = by the 
artificial 1/2 added to the angular momentum number 
I. Our wave function doesn't have the right power law 
behavior near the origin but there is only one divergence 
which is due to the right turning point. Hence, we see 
that the wave functions obtained from a systematic ex- 
pansion in powers of h without any ad hoc manipulation 
of the hydrogen problem are at least as accurate as those 
obtained from the problem with LM. 
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2. WKB wave functions and exact wave functions for 
;round state. The WKB wave functions diverge at the 
points, r is measured in units of the Bohr radius oq. 




FIG. 3. WKB wave functions and exact wave functions for 
the excited 3s and 3p states. Again the full, dashed, and 
dotted lines denote the exact, WKB(LM), and the WKB(SE) 
wave functions, respectively. 



Finally, we calculate radial dipole matrix elements be- 
tween states with angular momentum I and / ± 1 . Using 
the restricted interference approximation, we have 

Rin{E,l)^ (28) 
\§dr (r, E, L)r¥-^ (r, E + AE, L±n). 

Expanding this in powers of ?i, one finds for the leading 
order term 



Ri^\E^)^a, 



V? d 



An? de 



JA„(Ane) 



(29) 



± 



JA„(Ane) 



1/2 

where gq is the Bohr radius, e = [l — (^/n)^] the ec- 
centricity, and Jniz) a Bessel function. Naccache |l^] 
has obtained this leading order term from the Heisen- 
berg correspondence principle. The quantum correction 
of first order in h is found to read 



R 



±(1) 



iE,l) 



, 1±1 5 o±(0) 



(30) 



(EJ) 



1/2 

with the angular frequency llj{E) — [— 8£'^/(me^)] of 
the Kepler problem. In Tab. 1 the semiclassical dipole 
elements are compared with the exact ones for some spec- 
tral series. We note that for large n and I and small An 
the WKB results give rather accurate estimates of the 
exact values. This is expected from a semiclassical ap- 
proximation. 

In summary, we have shown that a systematic semiclas- 
sical expansion of the hydrogen problem about the Kepler 
problem yields remarkably accurate results. In contrast 
to the common belief no modification of the WKB ex- 
pansion is necessary when the centrifugal potential term 
is decomposed in the classical centrifugal potential and a 
quantum correction. The same method can be employed 
to other problems with radial symmetry. 
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n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Is-np 
2p-nd 
4p-ns 


1.090(1.290) 


0.512(0.517) 
4.542(4.748) 


0.339(0.306) 
1.816(1.71) 


0.257(0.209) 

1.104(0.975) 
4.673(4.600) 


0.208(0.155) 

0.802(0.662) 
1.864(1.788) 


0.177(0.121) 

0.641(0.492) 
1.120(1.044) 


0.154(0.098) 
0.543(0.386) 
0.794(0.718) 


0.137(0.082) 

0.478(0.314) 
0.614(0.539) 


0.124(0.069) 

0.432(0.263) 
0.501(0.427) 



TABLE I. WKB(SE) dipole matrix elements in units of Bohr's radius ao and exact quantum mechanical values in parenthesis. 
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